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We identify a universal mathematical structure in microscopic deterministic traffic models (with 
identical drivers), and thus show that all such existing models in the literature, including both 
the two-phase and three-phase models, can be understood as special cases of a master model by 
expansion around a set of well-defined ground states. This allows any two traffic models to be 
properly compared and identified. The three-phase models are characterised by the vanishing of 
leading orders of expansion within a certain density range, and as an example the popular intelligent 
driver models (IDM) is shown to be equivalent to a generalized optimal velocity (OV) model. We 
also explore the diverse solutions of the generalized OV model that can be important both for 
understanding human driving behaviours and algorithms for autonomous driverless vehicles. 
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It is of great interest, both theoretically and prac¬ 
tically, to understand via simple models the emergent 
behaviours of the complex systems containing a large 
number of interacting components. Examples like crowd 
dynamics pQ [5], highway traffic systems [31 |3], and the 
more recent urban traffic flows |S] have attracted physi¬ 
cists for decades. Despite the observed complex patterns 
from these systems, some of the essential characteristics 
are universally well-defined [6l[7|. In contrast to the tradi¬ 
tional many-body physical systems (e.g. involving iden¬ 
tical particles like electrons), the crowd or traffic systems 
lack almost any symmetry at the microscopic level: even 
individual components are different from one another, 
with intrinsic stochasticity. This poses great challenges in 
finding the simplest model with enough predictive power 
to adequately characterise these complex systems. 

The lack of symmetries at the microscopic level leads 
to a certain arbitrariness in the model construction. The 
modelling of the highway traffic system has led to a 
plethora of trafhc models Om inj|7] . Analysis of the em¬ 
pirical data has led to profound understanding of the traf¬ 
fic dynamics [51|S]; however due to non-linear interactions 
between drivers, simple interactions can lead to complex 
spatiotemporal patterns. It is thus difficult to decide 
which detailed driving behaviours or specihc functional 
forms should be used in the model. This is the primary 
concern behind the controversy between the two-phase 
and three-phase traffic theories [TOl [T3] . which leads to 
two classes of traffic models. While both the two-phase 
and the three-phase traffic models can simulate the scat¬ 
tering of the data points on the flow-density plane when 
the traffic is congested, the two classes are distinguished 
by the presence or absence of a continuous fundamen¬ 
tal diagram in the flow-density plotjT]. Moreover, the 
lack of a fundamental principle in understanding various 
different traffic models makes it difficult to decide if, in 
addition to the non-linear interactions between the vehi¬ 
cles, factors like stochasticity or diversity of driving be¬ 
haviours are also essential to certain observed empirical 


features. 

To tackle the dilemma of the apparent 
(over)abundance of traffic models, we need answers 
to the following questions: a). How do we properly 
characterise the differences between two traffic models? 
b). Is there a standard way of extending an existing 
traffic model or construction of new traffic models? c). 
Is there a standard way in selecting the best traffic 
model based on the empirical data? In this Letter, we 
answer in details the hrst two questions for deterministic 
microscopic models with identical drivers, by proposing 
a general framework of obtaining all such models from 
a master model. The controlled expansion around 
properly defined “ground states” of the master model 
shows that the two-phase and three-phase traffic models 
are both just special cases. The general framework also 
allows us to classify existing and new traffic model, and 
compare explicitly various types of approximations in a 
well-dehned way. 

The answer to the third question lies in the fact that in 
principle the master model can be obtained empirically, 
via a renormalization-like procedure by averaging over 
unimportant factors influencing the reaction of individual 
drivers. The details can be found in[T6] and will also be 
discussed elsewhere. We will now develop the general 
framework by assuming the following form of the master 
model 

~ f (1) 

Here the assumption is that the bumper-to-bumper dis¬ 
tance of the vehicle, hn, its velocity Vn, and the ap¬ 
proach velocity A.Vn = Vn+i — Vn are the important fac¬ 
tors affecting the acceleration of the n**' vehicle, and all 
other factors are averaged over. One can choose more 
(or different) dynamic variables in the function /, but 
for the purpose of illustration we choose this simple yet 
sufficient case. 

The natural time and length scale of the traffic system 
is given by the two statistically robust empirical quanti- 
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FIG. 1. Various possible sections of the acceleration (given by 
/ in Eq.|^) at different h, as a function of v and evaluated at 
Av = 0. The non-analyticity of c) can be well approximated 
by d). The specific form of such acceleration can in principle 
be obtained from empirical data|16| 


ties: pj, the density of the vehicles within a wide mov¬ 
ing jam, and Vjnaxj the maximum velocity (or the speed 
limit). Since we will carry out a controlled expansion 
of Eq.Q later, all quantities in Eq.Q are dimensionless 
after a chosen set of scales. While formally the model 
and its expansion are independent of the choice of the 
scales, for practical purposes one can choose either the 
length scale p~^, the time scale or other sets 

of the physically relevant scales, depending on different 
purposes and the particular system under study. 

The properties of the master model can be studied 
by plotting f {h,v, Av) as a function of v at fixed h in 
Fig.Q, where we take Av — 0 and focus around zero 
accelerations. The master model has a fundamental di¬ 
agram if for any value of h the acceleration decreases 
monotonically with v and crosses the x-axis only once, 
as shown in Fig.([^) or Fig.([^). This corresponds to 
the two-phase traffic models with a unique relationship 
between the flow and density. The difference between 
two-phase models as reflected by the detailed plots can 
be quantified by the proper Taylor expansions as we will 
show later. 

The plots in Fig. &) and Fig. 0) are cases where 
for certain range of h, there are more than one steady 
state (given by the vanishing acceleration) at the same 
density (^ h~^), but with different average velocities 
(given by the x-axis intersections). They correspond 
to the microscopic three-phase models proposed in the 
literature [m [T^. For Fig. the number of steady 
states at the same density equal to the number of inter¬ 
sections at the x-axis, though only the ones with negative 
gradient are stable. The Taylor expansions around differ¬ 
ent stable points where d = 0 leads to different optimal 


velocity functions within different velocity range. If for 
certain range of h the acceleration is not analytic in v as 
shown in Fig.([^), steady states can occur at the same 
average traffic density but for any velocities at which the 
acceleration is zero. Since the plot corresponds to the 
case where Av = 0, this is exactly the case of speed 
adaptation mechanism that leads to the scattering of the 
flow-density plot at the synchronized phase. Both these 
two cases are commonly used in various forms for the 
construction of three-phase models. 


The key assumption of the three-phase traffic theory is 
that the scattering of the flow-density plot in the “syn¬ 
chronized phase” corresponds to a multitude of steady 
states with non-unique flow-density relations [7]. it is 
however empirically difficult to verify if those states are 
indeed steady states. The macroscopic quantities like 
flow and density are empirically obtained from the aggre¬ 
gated raw data fmiT^ . A transient state with qualitative 
or quantitative features that evolves very slowly will be 
captured by the sensor as the “steady state”. Thus in 
most cases we cannot distinguish such long lasting quasi¬ 
steady states from the steady states. 


We thus argue that instead of dealing with a non- 
analytic function as shown in Fig.([^), we can for most 
practical purposes replace it with an analytic func¬ 
tion with an inflection point at zero acceleration (i.e. 
Fig.®)). Indeed the function in Fig. 0) can be approx¬ 
imated by an analytical function to any degree of accu¬ 
racy. Such a model still has a fundamental diagram, but 
for any state not too far away from the inflection point, it 
evolves very slowly with very small accelerations. Numer¬ 
ical simulation of such analytic models show very similar 
scattering of the flow-density plot and the spatiotemporal 
patterns as the proposed three-phase models within any 
physically reasonable time scale. Therefore fundamen¬ 
tally there is no strict boundary between the two-phase 
and three-phase models, and we will proceed to show the 
way to differentiate them is no different from the way of 
differentiating different two-phase models. 


We will now only deal with analytic master models. 
The ground states of the model are defined as the solu¬ 
tions with all vehicles equally spaced apart with headway 
h and travelling at the same time-independent velocity 
V. In Eq.Q the distance (or headway) is scaled by p~^ 
while the time is scaled by the tilded velocities 

and accelerations are derivatives of the scaled distance by 
the scaled time. Thus all quantities with tilde are scaled 
to be dimensionless. The ground states are indexed by 
the average density h~^, and at some densities there may 
be more than one ground states (i.e. Fig.®)). We can 
thus expand the master model around a chosen ground 
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state: 

Uji = ^ ^ ^p,q (j^n^ ^op (^) 

p,q 

ft "I _ 1 a-+’/ 

’’■■> y") ~ pk< arv„»iM,„ 
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where Vop (h) satisfies / {h, Vop (h) ,0) = 0. Clearly any 
existing two-phase models can be expanded this way, and 
they can be quantitatively compared by the set of coef¬ 
ficients Kp^q. The optimal velocity (OV) model[T7] is the 
special case where is a negative constant while all 
other Kp^q vanishes. The related FVD model has an ad¬ 
ditional constant ko,i[IH]i while the asymmetric FVD[TU] 
can be understood as having additional non-linear correc¬ 
tions with non-vanishing Ko,qj<l > 1- 

Within this framework the three-phase models and the 
phase transition can be understood as follows: there is a 
certain range of the density (typically when the bumper- 
to-bumper distance is smaller than some “synchroniza¬ 
tion gap” but larger than that within a wide moving jam), 
s:p,o vanishes for p < po so that the leading terms involv¬ 
ing Vn are of higher orders. The larger the value of poi 
the more long lasting the transient states are around the 
inflection points. Thus if one can obtain / empirically, 
the assumptions of the three-phase traffic theory can be 
tested microscopically. 

Seemingly different traffic models can be shown to be 
qualitatively equivalent under this scheme of controlled 
expansion. As an example we look at the IDM model 
which has the following form: 


(In — d 






h* {v, Av) = So + Si 



+ Tv + 


vAv 

2\fah 


( 4 ) 

( 5 ) 


Here we use dimensionful and physical quantities, with 
the standard set of parameters [l3j a — 0.73ms~^,b = 
1.67ms~^,vo = 33ms“^,so = 2m, si = 0m,T = 
1.6s, (5 = 4. The ground state is unique for any density; 
the relationship between h and Vop is given by 

= (so + KpT) h - M (6) 


The Taylor expansion around the ground state configu¬ 
ration has a finite number of terms: 

p—4,q—2 

an= ^ Ap,q {Vn - Vop {hn)Y Av^ (7) 

p=0,q=0 

In general the majority of the approach velocities of the 
highway traffic, as well as the deviation of the vehi¬ 
cle velocities from the optimal velocities, is on the or¬ 
der of lms“^, as verified from the numerical simulation 



FIG. 2. (Color Online.)The coefficients of expansion of the 
IDM in Eq.(^ as the function of the bumper-to-bumper dis¬ 
tance h, and Xpq is defined in Eq.Q and normalised to have 
the same dimensions 


of Eq.Q. We can thus de-dimensionalize velocities in 
Eq.0 in this unit, normalizing all the non-vanishing co¬ 
efficients of expansion to have the unit of the physical 
acceleration, and plot them in Fig.([^ as functions of the 
density. One can clearly see that the dependence of Xpq 
on traffic density is quite intuitive. The driver is only 
sensitive to Au„ at intermediate traffic densities. In the 
congested traffic the drivers are more alert (large Aio), 
and to keep a safe bumper-to-bumper distance is more 
important than to synchronize the velocity. 

It is thus probably enough from Fig.([^ to just keep 
AiO; Aoi, All for the qualitative features of the traffic dy¬ 
namics to be preserved. Numerical simulations show even 
ignoring An may be good enough. Such a truncated IDM 
model can also qualitatively simulate the approach dy¬ 
namics of a single vehicle rather well (see Figj^, even 
though in this special case the truncated terms are pre¬ 
sumably not small due to the large approach velocity and 
large deviation of the vehicle’s velocity from its optimal 
velocity. For the resulting truncated model, further tun¬ 
ing and optimization of the model are possible by fitting 
Aio and Aqi using simple functions (such as the step-like 
functions), to improve the quantitative agreements with 
the original IDM. 

Truncation of higher orders and simplification of the 
expansion coefficients can be applied to other models like 
Shamoto’s|20). We would also like to emphasize that it 
is not a priori true that the IDM or Shamato’s mod¬ 
els are more realistic than their simplified counterparts. 
The tuning of the expansion coefficients should be based 
on the experimental measurement of the master model. 
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FIG. 3. Here we simulate the dynamics of a single vehicle 
accelerating from zero velocity and approaching an obstacle 
located at 2.5km away from the starting point. The top two 
plots are the time evolution of the velocity and acceleration, 
simulated from the full IDM model. The bottom two plots are 
the time evolution of the velocity and acceleration, simulated 
from the linearlized IDM model by keeping only Aio and Aqi . 


and it is possible to achieve better agreement with the 
empirical observations even with the truncated models. 

The original OV models [T7H1T] are the simplest case of 
Eq. (§ where the coefficients of expansion are constants 
and only the term linear in (z;„ — Vop {hn)) is kept. The 
phase diagram of such models are well- known p2l42()j . To 
understand the more realistic traffic models as shown 
above, and for tuning and construction of new traffic 
models, it is also important to understand the emergent 
characteristics of the generalized OV models with general 
coefficients of expansion and higher orders of expansion. 
This is in general difficult due to the non-linear inter¬ 
action, and we will just explore some simple examples, 
leaving a more detailed discussion elsewhere. 

We start by generalizing the simplest artificial OV 
model as follows: 

a„ = A (/i„) (tanh (/i„) - w„) (8) 

It is useful to plot A together with the coexistence curve 
(CC) as well as the neutral stability (NS) line of the orig¬ 
inal OV model. The CC at /i < 0 is given by the mini¬ 
mum bumper-to-bumper distance of the cluster solution, 
which we call the /imin-branch; naturally the other half is 
called the /imax-branch[57]. The linear stability condition 
at average density is given by A (/iq) > 2sech^/io. 

The qualitative features of Eq. ® can be roughly clas¬ 
sified by the way A intersects witn the CC and the NS 
line. If A intersect with the /iminr and the /imax-branch at 
most once separately (see Fig^^|^)), the constant density 
{hn = const) solution of Eq.® can be stable or unstable 
against the development of cluster solutions with mini¬ 
mum and maximum bumper-to-bumper distances /imin 
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FIG. 4. (Golor Online.)A11 the insets are plots of headways as 
a function of the car index, a). A^^\ A*-^^ and A^^l are examples 
of different functions of A in Eq.Q. For cluster solutions the 
maximum or minimum space gaps can be uniform as c)., from 
aC) and A^^>), or fluctuating as d)., from A^^l). b). A^^^ and 
A(2) 

are examples of different functions of A. The vehicle 
density corresponding to the thickened part of A, if unstable, 
leads to regular cluster solutions, e). shows incomplete cluster 
solutions when the vehicle density satisfies A (/iq) < 2sech/i0i 
which is linearly unstable. 


and /imax, qualitatively the same as the original OV 
model; both /imin and /imax are independent of tlw den¬ 
sity of the traffic. Clearly only A S [A (/imin) , A (/imax)] 
contributes to the cluster structure, with the following 
relationship: 

^max/min ^ [^max/min (A {h 

H-ii n. ^ ^ 1 rxitL3c ! imi H {\{h 

max ))] (9) 

If A (/imin) ^ 2sech(/imin) A (/i-max) ^ 2sech(/lmax): the 
vehicles in the clusters (or anti-clusters) are also linearly 
unstable and their burnper-to-bumper distances will fluc¬ 
tuate (see A^^^ in Fig.(|4^) and the inset), in contrast to 
the original OV model. 

If A intersects with the /imin-branch (/imax-branch) 
more than once (Fig.(|^)), the traffic dynamics can be 
different at different vehicle density. For /iq within the 
thickened part of A*-^^ and A^^\ the stability of the 
hn = const solution and the development of the cluster 
solution are similar as the cases in Fig.([^). For other 
part of ho, if A(/io) < 2sech(/io), the /i„ = const solu¬ 
tion is linearly unstable but clearly the cluster solution 
cannot be fully developed due to the conservation of the 
vehicle density. Partially developed cluster solutions with 
strong oscillations are obtained as shown in the inset of 
Fig. I^) EH. 

Higher orders of {Vop — Vn) generally shift the values 
of /imin and /imax- Here we look at simple cases with 
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FIG. 5. a). Cluster solutions in the unstable phase, when 
higher order terms in {v — Vop) are just small perturbations 
(Top left), b). Multiple equilibrium densities at the same 
velocity, corresponding to the case in Fig.([^), when the 
h„ = const solution at each x-axis intersection with negative 
gradient is stable (Top right), c). The unstable phase when 
only (u — Vopf' term is present in the expansion with constant 
expansion coefficient (Bottom left), d). The unstable phase 
when only {v — Vop)^ term is present in the expansion with 
constant expansion coefficient (Bottom right). 


constant coefficients of expansion. For situations like the 
ones in Fig. the stable phase can have multiple equi¬ 
librium densities with vehicles all traveling at the same 
velocity (see Fig.([^)). For three-phase models with van¬ 
ishing low order terms (at least the linear oder), the NS 
line is no longer applicable. Since the accelerations are 
very small even when the vehicles deviates from their 
equilibrium density, it is very difficult for the clusters to 
form, as shown in Fig.([^) and Fig.([^). 

One should note that our method can be applied to 
more general cases, in particular with more realistic op¬ 
timal velocity functions. The cluster solutions, the CC 
and NS line are universal features for the more compli¬ 
cated models. Though step functions of A are used in 
many of the arguments above for its simplicity, A can be 
deformed moderately or smoothened without altering the 
qualitative features. Since only A S [A (^min) ) A (^max)] 
determines the instability of the model, one can alter 
the shape of A outside of that range arbitrarily without 
changing the cluster structure or the phase diagram. For 
example, the inset of Fig.(|^) will be realized by a mono- 
tonically decreasing A intersecting the /i„iax-branch of the 
CC twice, as is the case for the IDM model. 

In conclusion, we propose that microscopic traffic mod¬ 
els should start with identifying the proper optimal ve¬ 
locity function, from which the family of ground states 
are defined. Simple traffic models are constructed as spe¬ 


cial cases of the expansion around those ground states, 
with coefficients of expansion in general dependent on the 
vehicle density. In particular the “synchronized phase” 
can be formally understood as the result of the mod¬ 
els where the leading order expansion about the optimal 
velocity vanishes within some density range. The theo¬ 
retical framework also proposes standard ways of model 
construction, and validation by the empirical data. The 
understanding of various possible solutions to a more gen¬ 
eral class of OV models could also be useful for designing 
the algorithms of autonomous driverless vehicles for large 
scale highway transportation. 

We thank Prof. Dirk Helbing for useful discus¬ 
sions. This research was partially supported by Singa¬ 
pore A*STAR SERC “Complex Systems” Research Pro¬ 
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